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Abstract 

The cross sections of the processes 4 He(7,p) 3 H and 4 He(7,n) 3 He are calculated taking into 
account the full final state interaction via the Lorentz integral transform (LIT) method. This is 
the first consistent microscopic calculation beyond the three-body breakup threshold. The results 



are obtained with a semirealistic central NN potential including also the Coulomb force. The cross 



sections show a pronounced dipole peak at 27 MeV which lies within the rather broad experimental 
band. At higher energies, where experimental uncertainties are considerably smaller, one finds a 
good agreement between theory and experiment. The calculated sum of three- and four-body 
photodisintegration cross sections is also listed and is in fair agreement with the data. 

PACS numbers: 25.20.Dc, 21.45.+V, 24.30. Cz, 27.10. +h 
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I. INTRODUCTION 



The photodisintegration of the 4 He nucleus has a longstanding history. First experiments 
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on the (7, n) reaction were performed some 50 years ago |l|. In the following 25 years various 
experimental groups carried out measurements for both the 4 He''" - ' 
reaction channels and the inverse capture reactions 
Dramatically conflicting results have been obtained as a result of this work. The (7,p) data 
were consistent in showing a rather pronounced resonant peak close to the three-body 
breakup threshold. At the same time, the (7, n) data at low energy were very spread, 
and measurements showed either a strongly pronounced or a rather suppressed giant dipole 
peak. In 1983 a careful and balanced review of all the available experimental data for the 
two mirror reactions was provided Q|. A strongly peaked cross section at low energy for 
the (7,p) channel and a flatter shape for the (7,n) one was recommended by the authors. 
Three new experiments on the (7,p) reaction were subsequently carried out, two of them 
contradicting and one confirming the recommended cross section. In particular in 131 and 

n n 

|18| a suppressed dipole resonance was found, whereas [19J proved to be in agreement with 
the previous strongly peaked results. Measurements of the ratio of the (7,2?) to (7, n) cross 
sections in the giant resonance region were performed as well and, at variance with the cross 



Finally, 



sections recommended in Q], results very close to unity were reported [21 
in 1992 additional cross section data were deduced from a Compton scattering experiment 
on 4 He jj^o]]. A strongly peaked cross section for the 4 He total photoabsorption was found 
suggesting a (7,^) cross section considerably larger than the one recommended in jl6|. 

In 1996 the first theoretical calculation of the two-fragment breakup cross section with 
inclusion of final state interaction (FSI) was performed in the energy ran ge b elow the three- 



fragment breakup threshold 231 ] . The semi-realistic MTI-III potential 24] was employed. 
The results of [23] showed a rather suppressed giant dipole peak. The following year a 
calculation of the total 4 He photoabsorption cross section up to the pion threshold was 
carried out . Full FSI was taken into account in the whole energy range via the Lorentz 
integral transform (LIT) method (2^ that does not require calculating continuum wave 
functions. The four-nucleon dynamics was described with the same NN potential as in j^J . 



Different from the previous work a pronounced giant dipole peak was found 25]. These 



results have been reexamined in |2j|. A small shift of the peak position has been obtained, 



2 



but the pronounced peak has been confirmed. 

In order to understand the origin of the large total photonuclear cross section in terms 
of the various channel contributions, a separate calculation of the channel cross sections 
is necessary. The present paper reports a theoretical study of the two-fragment breakup 
processes 4 He(7,p) 3 H and 4 He(7,n) 3 He. For these exclusive reactions we investigate the 
issue of the giant dipole peak height, and we also calculate the cross section at higher 



energies, up to the pion threshold. The calculation inc 
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udes full FSI via the LIT method 



[28] and employs the semi-realistic MTI-III potential 

The LIT method has already been successfully tested for exclusive reactions in the 
d(e,e'p)n process j^. The present calculation represents the first application to A > 2. 
We would like to emphasize that in our study FSI is taken into account rigorously also be- 
yond the three-and four-body breakup thresholds. In addition, combining our cross sections 
with the total photoabsorption cross section calculated in 2Jj,|27[ for the same potential, we 
obtain the sum of three- and four-body photodisintegration cross sections of 4 He, for which 



some experimental data are available 



II. GENERAL FORMALISM 



A. Cross Section 



The total exclusive cross section of the 4 He photodisintegration into the two fragments N 
and 3, where iV refers to the scattered proton (neutron) and Prefers to the 3 H ( 3 He) nucleus, 
is given by 

°n,3 = Yc^&T j \(^n, 3 ( En ^ \ D *\ ^ a )\ 2 dtt k , E N)3 = to^ + E a . (1) 
In the equation above we neglect the very small nuclear recoil energy. With fi and k we 
denote the reduced mass and the relative momentum of the fragments, respectively, co> 7 is 
the incident photon energy, ^ a is the ground-state wave function, ^Jj 3 is the final-state 
continuum wave function of the minus type pertaining to the N,3 channel and E n anc 
i?7V 3 are the energies of the corresponding initial and final states, respectively. As in |25l 127 
only transitions induced by the unretarded dipole operator, 
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^ = E r ^- ( 2 ) 

3=1 



are taken into account. Here rj denotes the third component of the j-th nucleon isospin 
and Zj represents the z component of the distance between the j-th nucleon and the center 
of mass of the system. This form of the transition operator already includes the leading 



contributions to the 
35lQ). 



effects of meson exchange currents via Siegert's theorem. Additiona 
total cross section are small at the energies considered here (see also 

The main difficulty in the calculation of such a cross section is the presence of the con- 
tinuum wave function ^^3(^^,3) i n the transition matrix element 

T N , 3 (E N , 3 ) = (^3(^,3) \D Z \ *«) ■ (3) 

As it will be explained in the next subsection, with the LIT method one is able to perform 
an ab initio calculation of this transition matrix element in the whole energy range below the 
pion threshold without dealing with the continuum solutions of the four-body Schrodinger 
equation. 

B. The LIT Method for Exclusive Reactions 

For the case of an exclusive perturbation-induced process, all the information about 
the reaction dynamics is contained in the transition matrix element of the perturbation O 
between the initial (^/q) and final fifj) states, 



Tf(E f ) = (*j(E f ) O • (4) 
The calculation of such a matrix element can be carried out with the LIT method as outlined 



in the following [28|, |29( . For simplicity we restrict our discussion to an exclusive reaction 
leading to a final state with two fragments, a fragment a with n a nucleons and a fragment 
b with rib = A — n a nucleons, though more complex fragmentations of the initial A-body 
system can be treated as well. Denoting with H the full nuclear Hamiltonian we have a 



formal expression for ^ f =a ^ in terms of the "channel state" 4 > f =a b{Ef= a ,b) 
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where A is an antisymmetrization operator. When there is no Coulomb interaction in the 
system, the channel state (p~ b (E a ^) is the product of the internal wave functions of the 
fragments and of their relative free motion. Correspondingly, V in Eq. (JBJ) is the sum 



of all interactions between particles belonging to different fragments. In the general case, 
(f)~ h {E a ^) is chosen to account for the average Coulomb interaction between the fragments, 
and the plane wave describing their relative motion is replaced by the Coulomb function of 
the minus type. We write down this function in the partial wave expansion form 



$ a (l,...,w a )$ h (w a + 1,...,A) 
(2tt) 3 / 2 



i' e -«i(*) 



we(k;r) 
kr 



Y im (n r )Y; m (n k ). (6) 
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Here $ a (l, n a ) and $b(n a + 1, A) are the internal wave functions of the the fragments, 
r = (r, fi r ) = R" m — R^ m represents the distance between them, and the energy of the 
relative motion is /c 2 /2/z = E a ^ — E a — Ef,, where E a and E b are the fragment ground state 
energies. The functions we(k;r) are the regular Coulomb wave functions of order £, and 
5e{k) are the Coulomb phase shifts j^. 

The internal wave functions of the fragments are assumed to be antisymmetrized and 
normalized to unity, so that the properly normalized continuum wave function in Eq. (jSJ) is 
obtained via application of the antisymmetrization operator 



A 



n a \n b \ 
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(7) 



E E r« 
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where Vij are particle permutation operators [34]. 

When one inserts Eq. (jSJ) into Eq. (j3J) the transition matrix element becomes the sum of 
two pieces, a Born term, 



a,b {Ea,b) 



KM(Ea,l 



and an FSI dependent term, 
T a T(E a , b ) 



J a,b 



(Ea,b) 



VA 



A O 



1 



E a h + is — H 



;0 



(8) 



(9) 



While the Born term is rather simple to deal with, the determination of the FSI dependent 
matrix element is rather complicated. We treat this term within the LIT approach as 
outlined in the following. 

Let ty u (E u ) be the eigenstates of the Hamiltonian labelled by channel quantum numbers 
v and normalized as ^ v \^ v i) = 5(v—u'). Using the completeness relation of the set ^> V {E V ) 
the matrix element T^ b SI (E ai b) can be written as: 

1 



rpFSI 

1 a,b 



(Ea,l 



du{<j>- b (E aib )\VA\* v (E„)) 
3 F a ,b(E) 



e~, E a ,b + ie — E 



-dE 



E a ,b + ie — E v 
- tTTF a ,b(E a , b ) + V 



* u (E v )\0\*o) 

Fa,b(E) 



dE , (10) 



where F a>b (E) is defined as 

F nh {E) = yfjvL,{E a , b ) 



VA 



<S> V (E V )){^ V (E V ) O * )6{E-E, 



and E t h is the lowest excitation energy in the system i.e. the breakup threshold energy. 
To calculate T a ^ one needs to know the function F aib for all energy values. The direct 
calculation of F a ^ is of course far too difficult, since one should know all the eigenstates ty u 
for the whole eigenvalue spectrum of H. However, an indirect calculation of F a ^ is possible 
applying the LIT method. To this end, one introduces an integral transform of F a ^ with a 
kernel of Lorentzian shape, 

F a>b (E) JB1 [ x F a , b (E) 



L [F a ,b] (a) 



{E _ aR) 2 + a 2 



dE 



^du(cf>~ b (E a>b 

V 2 (o) I $i(c7 



VA 
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*„(£„)) U V (E V 



dE 
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with a = a R + iaj. Here we denote 



^(a) = (H — a^O |* > , ^2(0-) = (H- a^AV^AE^)) . 



(12) 



(13) 



Equation (JT2*Jl shows that L [F a>b \ (a) can be calculated without explicit knowledge of F ajb 
provided that one solves the two equations, 



(H-a) 
(H-a) 



Vi) = O|* ) , 
9 2 ) = AV\K b (E a , b )) 



(14) 
(15) 



which differ in the source terms only. 

It can be seen that ^ and ^2 are localized i.e. have finite norms. Assuming that 
07 7^ this is ensured by the fact that the source terms in the right-hand sides of Eqs. 
()14j) and ()15|) are localized. In fact in the coordinate representation the source term in Eq. 
(JUjl decreases exponentially for increasing distances between particles. The source term 
in Eq. (|15|) contains two potential contributions i.e. the nuclear force component and the 
Coulomb interaction. The contribution due to the nuclear interaction vanishes exponentially 
for increasing distances between particles. Since the average Coulomb potential between the 
fragments is already taken into account via Eq. ©, the remaining Coulomb contribution 
behaves like r~ 2 , as the distance r from the fragments increases. Therefore the source term 
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remains localized also in the presence of the Coulomb force. In our calculation this Coulomb 
term arises in case of the p, 3 H channel and its contribution to the results is found to be very 
small. 

In case of breakups into more than two fragments it is convenient to write Eq. (|12|) in 
another form: 



L [FaA {<?) 



1 



2iai 
1 

2iai . 



VA 



H-a H-a* 



O 



- b (E ajb )\VA ~ (fc, b (Ea,b)\VA 



(16) 



Here ^i(cr) is the same as above, and ^[(a) satisfies 

(ff-O |*1(*)) = 6 1* > . (17) 

Therefore, also in this case, one has to deal with a localized solution. 

When solving Eqs. (fTH) . (fTKj) and (fTTj) it is sufficient to require that the solutions are 
localized, and no other boundary conditions are to be imposed. Therefore, similar to what is 
done in bound state calculations, one can use an expansion over a basis set of localized func- 
tions. A convenient choice is the basis set consisting of correlated hyperspherical harmonics 
(CHH) multiplied byhyperradial functions. 

As discussed in [27] for the case of the total 4 He photoabsorption cross section, special 
attention has to be paid to the convergence of such expansions. A rather large number of 
basis states is necessary in order to reach convergence, thus leading to large Hamiltonian 
matrices. Instead of using a time consuming inversion method we directly evaluate the scalar 
products in |T6|) with the Lanczos technique [37]. Inserting the Lanczos orthonormal basis 
{(fi,i = 0, ...,n} into Eq. (|TH|) . where (p a is taken to be the right-hand side (normalized to 
unity) of Eqs. (JH|), (fTTjl . one gets 



L[F atb ]{a) 
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H 
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(18) 



The matrix elements (ifi\(H — a)~ l \(p ) can be written as continued fractions of the Lanczos 
coefficients. 

After having calculated L[F a)b ](a) one obtains the function F atb (E), and thus T ajb (E a>b ), 



via the inversion of the LIT, as described in 



38|. 
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III. RESULTS 
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The ground states of 4 He, 3 He and 3 H as well as the LIT in Eq. (|18|) are calculated 
using the CHH expansion method. In order to speed up the convergence, state independent 
correlations are introduced as in 3^|. We use the MTI-III j^J potential and identical CHH 
br the ground state wave functions of 4 He and of the three-nucleon systems as 
respectively. 

Because of the choice of the excitation operator (see Eq. (J2J), in solving Eqs. (jHj) and 
()15|) the HH must be characterized by the quantum numbers L=l, S=0 and T=l. In this 
calculation the maximal value of the grand-angular quantum number, K max , is equal to 29. 
Such a high value of K max has been made possible by neglecting states which have proved 
to give very small contributions to the LIT. The choice of these states has been done in the 
following way. Our hyperspherical harmonics (HH) depend on three Jacobi vectors £i, £2, £3, 
among which £1 represents the distance between a pair of particles. The HH entering the 
calculation are obtained via symmetrization of the HH possessing a definite relative orbital 
momentum l\ associated with £1 and a definite grand-angular momentum K in the subspace 
of vectors £2, £3- We have found that beyond K max = 13 we can neglect all symmetrized 
HH which originate from the non-symmetrized HH with i\ greater than 2 and K greater 
than 3 and, beyond K max = 19, those with i\ greater than and K greater than 1. This 
selection is similar to that justified and used in non-correlated HH bound state calculations 
(see 3| and references therein). An additional selection has been performed with respect to 
the permutational symmetry types of the HH. For S — 0, T = 1 quantum numbers, HH of 
two permutational symmetry types enter the expansion, those belonging to the irreducib 



representations [/] = [+] and [/] = [— ] of the four-particle permutation group S4 |39j,|41|. 
For basis states constructed from the HH of the [+] type the probability for a nucleon pair to 
be in a relative even state and in particular in the s-state is substantially higher. We have 
found that the contribution to the LIT of states with the spatial symmetry [— ] is suppressed, 
and the corresponding HH with K values higher than 9 can be neglected in the calculation. 

In Fig. [H the convergence of the LIT with respect to K max is shown. One sees that for 
K max = 29 a good convergence is reached. 

In Fig. |21 we present our results for the 4 He(7,n) 3 He cross section together with exper- 
imental data. The difference between the total cross section and its Born approximation 
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shows large effects of FSI. At small energies FSI enhances the transition matrix element since 
the values of the final state continuum wave function inside the reaction zone become larger 
due to the attraction between the fragments. At higher energies the transition matrix ele- 
ment including FSI becomes smaller than the Born one. This is due to shorter wavelengths 
(higher particle momenta) leading to faster oscillations in the integrands. 

As it was pointed out in the introduction the experimental results do not show a unique 
picture. In the e nerg y region beyond 35 MeV there is an overall agreement of the data of 
4], ^SEll^l^BL whereas in the dipole resonance region big discrepancies are present. 
The data from (4 , UgI Q] show a rather pronounced dipole peak, while the more recent ones 
from show a flatter behavior (for more detailed information see also reviews of (7, n) 



experiments 



12 



peaked data of 



16 



). In the low-energy region our full calculation favours the strongly 
4J, |8 , 0, Q] . It is seen from the figure that the higher energy experimental 
cross section agrees quite well with our full calculation. We would like to point out the large 
effect of FSI in the high-energy tail. This is probably the region where large dependences 
of the cross sections from the details of the force and in particular of the three-body force 
are expected (see and where it would be desirable to have more accurate data. 

The present 4 He(7, n) 3 He results are at variance with the previous calculation carried out 
up to the three-body breakup threshold (AGS calculation with a rank 1 separable potential 
representation of the MTI-III potential |23||). 

The 4 He(7, p) 3 H cross section is compared with data in Fig. El The experimental situation 
is similar to the (7, n) case considered above. Again some experimental results show a 
pronounced giant dipole peak at low e nerg ies, as in case of , and 

others present much less strength, like |lj and jl8j |. At higher energies there is quite a 
satisfactory agreement among experimental data of different groups. As in the 4 He(7, n) 3 He 
case, our full results have a large cross section in the dipole resonance region and favour the 
strongly peaked data, while beyond the peak one finds a rather good agreement between 
our calculation and all data. 

In Fig. 13] we show the effect of the Coulomb interaction in the 4 He(7, iV)5 reactions. 
Though the effect is generally rather small, one notices a suppression of the cross section 
near threshold due to the Coulomb interaction in the (7,^0) case, compensated by an increase 
at higher energies. 

Since the 4 He(7, d)d reaction is not induced by the dipole operator of (j2J, the sum of 
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the 4 He(7,p) 3 H and 4 He(7,n) 3 He cross sections has to be equal to the 4 He total inclusive 
photoabsorption cross section below the n + p + d breakup threshold. Comparison of the 
sum of our exclusive cross sections in this region with the total inclusive photoabsorption 
cross section calculated independently can therefore serve as a test of our results. To this 
aim, we have calculated also the total photoabsorption cross section. This was done using 
the fact that the norm (^i(a) \ ^i(cr)\ represents the LIT of this total cross section. Our 
present CHH calculation reproduces the results obtained for the total photoabsorption cross 
section with the effective interaction HH method 27]. This can be considered as one more 
test of our calculation, and this is accomplished thanks to the present efforts made in order 
to increase the K max value up to 29. 

The comparison of the sum of our cross sections for the 4 He(7,p) 3 H and 4 He(7, n) 3 He 
reactions with the CHH total photoabsorption cross section below the three-body breakup 
threshold is presented in Fig. |S| The agreement of the two curves in this region is quite 
satisfactory. There is only some discrepancy very close to threshold. The discrepancy could 
presumably be resolved by a calculation of the LIT in the threshold region with a smaller 
(Tj. However, this also requires an increase of the number of hyperradial basis functions 
along with a more precise determination of the nine-dimensional integrals (at present Monte 
Carlo integrations), and thus a considerable improvement of the numerical calculation. This 
uncertainty is in any case smaller than the error bars of available experiments in that energy 
range. 

Figure also shows an indirect determination of the 3,4-body breakup cross section, 
which is obtained by subtraction of the 4 He(7,p) 3 H and 4 He(7,n) 3 He cross sections from 
the total photoabsorption cross section in the region where significant differences are found. 
This more-body breakup cross section is shown in Fig. |U] together with the three- and 
four-body disintegration data. Though there are some differences at lower energies and the 
peak of the theoretical cross section is more pronounced, one finds an overall fair agreement 
between theory and experiment. 

In conclusion we summarize our work. We have calculated the total photodisintegration 
cross section of 4 He in the two mirror channels p, 3 H and n, 3 He with the MTI-III potential 
using the LIT method and the CHH expansion. For the first time a microscopic calculation 
for these reactions has been carried out taking fully into account final state interaction also 
beyond the three- and four-body breakup thresholds. The pronounced dipole resonance 
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found in the total photonuclear cross sections 



25 



271 ] is almost exhausted by the two-body 



break up channel. The 4 He(7,p) 3 H and 4 He(7,n) 3 He cross sections we obtain have a very 
similar structure. Both show a pronounced giant dipole resonance peak at low energy as 
obtained in a previous calculation of the total 4 He photoabsorption cross section with the 
same NN potential model. This low-energy behavior of the cross section is in agreement 
with the ex— data offlflflfiflflHQ whereas other _ent S 
[12I Ibsl . I17I Ibst show less strength. At higher energies there are much less differences among 
the experimental data and one finds a rather good agreement with the theoretical results. 
Our indirect determination of the sum of three- and four-body breakup cross sections shows 
an overall fair agreement with the sum of three- and four-body break-up data. 

We hope to have convinced the reader that further theoretical and experimental efforts are 
necessary for a better understanding of the 4 He photodisintegration. The present theoretical 
results are fully consistent and complete regarding the treatment of the dynamics in the 
initial and final state, though obtained with a simple NN potential. In order to confirm 
the size of the giant dipole peak one would need a more realistic description of this process 
considering modern realistic NN potentials together with a three-nucleon force. 
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FIG. 2: 4 He(7, n) 3 He cross section up to 120 MeV (a) and 35 MeV (b): full result with FSI included 
(solid curve), Born approximation only (dashed curve); experimental data from 4] (full squares) 
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FIG. 3: 4 He(7,p) 3 H cross section up to 120 MeV (a) and 35 MeV (b): full result with FSI included 
(solid curve); ex peri mental data from [J] (full squares), [j| (full diamonds), [(| (open diamonds), j^| 

Ben squares), |ll| (downward open triangle s) , \l\ (open circles), Q| (upward open triangles), 
(upward full triangles), Q| (full circles), Ql (stars) and Q] (downward full triangles). 
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